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Abstract: Site measurements were collected at Mount John University Observatory in 2005 and 2007 
using a purpose-built scintillation detection and ranging system. Cf^{h) profiling indicates a weak layer 
located at 12 - 14 km above sea level and strong low altitude turbulence extending up to 5 km. During 
calm weather conditions, an additional layer was detected at 6 - 8 km above sea level. V{h) profiling 
suggests that tropopause layer velocities are nominally 12-30 ms~^, and near-ground velocities range 
between 2-20 ms~^, dependent on weather. Little seasonal variation was detected in either C^{h) 
and V{h) profiles. The average coherence length, ro, was found to be 7 =b 1 cm for the full profile at a 
wavelength of 589 nm. The average isoplanatic angle, ^o, was 1.0 ±0.1 arcsec. The mean turbulence 
altitude, /lo, was found to be 2.0 ±0.7 km above sea level. No average in the Greenwood frequency, /g, 
could be established due to the gaps present in the V(h) profiles obtained. A modified Hufnagel- Valley 
model was developed to describe the C%{h) profiles at Mount John, which estimates ro at 6 cm and 
Oo at 0.9 arcsec. A series of V{h) models were developed, based on the Greenwood wind model with 
an additional peak located at low altitudes. Using the C%{h) model and the suggested V{h) model for 
moderate ground wind speeds, fa is estimated at 79 Hz. 
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1 Introduction 

Astronomical images taken by ground based telescopes 
are subject to distortion caused by atmospheric turbu- 
lence. Adaptive optics (AO) provides a real-time solu- 
tion to compensate for an aberrated wavefront through 
the use of deformable optics in a closed-loop system. 
Accurate measurements and models of atmospheric tur- 
bulence are an essential tool in the design and optimi- 
sation of an AO system (Avila et al. 2001). 

Key parameters for the design of an AO system 
include the turbulence coherence length, ro, the iso- 
planatic angle, ^o, and the Greenwood frequency, fo- 
ro describes the effective telescope diameter for which 
nearly diffraction-limited resolution can be obtained 
if no attempt is made to compensate for atmospheric 
turbulence. It is defined as (Tyson 1991) 



0.423A;^ sec(C) / C%{h)dh 
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where /c = 27r/A is the wavenumber for a given wave- 
length A and C is the zenith angle. C%{h) is the re- 
fractive index structure constant and is a measure of 
the strength of a turbulent layer located at altitude h. 

The isoplanatic angle, Oo, describes the maximum 
angular separation between two objects for which tur- 
bulence induced distortions are essentially identical and 
is defined as (Parenti and Sasiela 1994) 



Oo 
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Unlike ro, ^o is dependent on h^^^ indicating that weak 
high altitude layers have a significant impact on ^o- 

Atmospheric turbulence is in a constant state of 
motion. The Greenwood frequency, /g, describes the 
rate at which the turbulence structure above a site 
changes with time. It is defined as (Tyson & Frazier 
2004) 



fo = 0.255 



k secC 
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where V{h) is the average horizontal wind velocity as 
a function of altitude h. fo determines how quickly an 
AO system must respond to adequately compensate for 
the aberrations induced by atmospheric turbulence. 

The 1-m McLellan telescope at Mount John Uni- 
versity Observatory (MJUO), located at Tekapo, New 
Zealand, is used for a variety of different astronomi- 
cal research and is known to regularly experience poor 
seeing (> 2 arcsec angular resolution) by observers 
(A. Gilmore (MJUO) 2006, private communication). 
This work is part of a feasibility study on installing an 
AO system to improve photometric images with the 
CCD photometer head and to improve light through- 
put into the HERCULES echelle spectrograph (Hearn- 
shaw et al. 2002) currently installed on the 1-m tele- 
scope. The elevation of MJUO is 1024 m. 

This paper discusses the C%{h) and V{h) profile 
measurements taken at MJUO and the models devel- 
oped for use in an AO system design. Section 2 outlines 
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the techniques used to measure C%{h) and V{h) pro- 
files and the purpose-built system. Sections 3 and 4 
discuss the data collected and the trends noted in the 
profiles obtained. Section 5 introduces the models de- 
veloped to describe atmospheric turbulence at MJUO. 
Concluding remarks are in section 6. 

2 Measuring Turbulence 

SCIDAR (scintillation Detection And Ranging) is a 
remote sensing technique that has been used at many 
different sites around the world to characterise opti- 
cal turbulence (Avila et al. 2001; Garcia-Lorenzo et al. 
2009; Masciadri et al. 2010; Prieur et al. 2001; Tokovinin 
et al. 2005). It uses the spatio-temporal covariance 
functions obtained from a sequence of short exposure 
images of the scintillation pattern seen at a telescope 
pupil to infer the C%{h) and V{h) profiles present 
above a site (Kliickers et al. 1998). 

SCIDAR measurements are commonly taken using 
a double star system, as indicated in Figure 1. Light 
from each star passes through the same region of a tur- 
bulent layer forming identical, but separated, scintil- 
lation patterns. The distance between the two scintil- 
lation patterns is directly proportional to the angular 
separation of the double star system, 0, and the height 
of the turbulent layer above the measurement plane, 
hi. 

Pupil-plane SCIDAR measures scintillation 
patterns seen at the telescope aperture. In doing this, 
a clear picture of optical turbulence in the free at- 
mosphere can be obtained. As scintillation is propor- 
tional to h^^^ (Roddier 1981), any scintillation result- 
ing from near-ground turbulence (NGT) is not readily 
detectable. Using a simple lens change, the measure- 
ment plane can be shifted to a virtual plane located at 
d below the telescope. U Hl is the height of the layer 
above the telescope then the height of the layer above 
the measurement plane becomes hi = — d\, where 
d is negative due to sign conventions. This increased 
propagation distance allows for scintillation from NGT 
to be adequately measured. This version of SCIDAR is 
known as generalised SCIDAR (Kliickers et al. 1998). 

The University of Canterbury SCIDAR system 
(UC-SCIDAR) is a purpose-built instrument designed 
to measure C]^(h) and V(h) profiles (Johnston et al. 
2004, 2005; Mohr at al. 2006, 2008a, 2008b). UC- 
SCIDAR saw first light on the McLellan 1-m tele- 
scope at MJUO in late 2003, at an approximate cost 
of USD$4000. To keep costs low the system design 
utilised primarily off-the-shelf components. Over the 
years, the system has evolved through several itera- 
tions. 

The current system consists of two channels, each 
with its own CCD camera and field lens (Figure 2). 
Light from the telescope is split using a 50/50 inten- 
sity beamsplitter (BS) to minimize the differences be- 
tween the two channels. The straight path, typically 
used for pupil-plane SCIDAR, consists of a fl2.7mm 
achromat lens (LI) mounted in a lens tube a focal 
length away from a Dragonfly Express CCD camera 
from Point Grey Research (CI). This camera uses a 




Figure 1: The concept of double star SCIDAR. 
Light from each star passes through the same tur- 
bulent region forming identical scintillation pat- 
terns separated by a distance proportional to the 
double star separation (j) and the height of the tur- 
bulent layer above the measurement plane hi = 
\hL — d\. Due to sign conventions d is negative. 




C1 



(a) Physical layout (b) Optical layout 

Figure 2: (a) Physical and (b) optical layout of 
the UC-SCIDAR instrument. See Section 2 for a 
detailed description of the system. 



Kodak KAL0340DM sensor which has a 640 x 480 grid 
of 7.4 /im square pixels. UC-SCIDAR has been config- 
ured to capture images at full resolution with a frame 
rate of 60 Hz. Operational exposure times range from 
0.5 to 5 ms. The lens used provides a nominal spatial 
sampling, Ar, of 1/125 m.pix"^ when the 1-m tele- 
scope is operating at a focal ratio of F/13.5. 

A second identical CCD camera is mounted in the 
side path (C2). As this channel is typically used for 
generalised SCIDAR, a flOmm achromat lens (L2) is 
mounted a focal length of LI away from C2. This pro- 
vides a measurement plane at approximately 3.5 km 
below the telescope. Figure 3 shows a typical pupil- 
plane and generalised SCIDAR scintillation image ob- 
tained using UC-SCIDAR. 

Assuming that exposure times and frame rates are 
short enough that turbulent elements move without 
distortion, then the 2D spatio-temporal covariance 
function for a double star can be written as (Avila 
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(a) Pupil-plane frame (b) Generalised frame 

Figure 3: Typical scintillation images from (a) 
pupil-plane and (b) generalised SCIDAR. 



et al. 2001) 

n 

Cb{p, 0, At) = J2iaCs{p - V{h^)At) 

+hCs{p-V{hi)M-(j)h^) 
+hCs{p - V{hi)At + 0/1,)}, (4) 

where is the angular separation of the double star, 
At is the time delay between consecutive scintillation 
images, n is the number of discrete turbulent layers 
present and Cs{p — V{hi)At) is the spatio-temporal 
covariance of a single star for the radial coordinate, p, 
in the direction of the double star. The coefficients a 
and b are given by (Avila et al. 2001) 



l + g^ 
(1 + a) 



,a = 10" 



(5) 



where Am is the magnitude difference between the 
double star components. Cb{p, (j), At) describes a se- 
ries of triplets where each corresponds to a different 
turbulent layer. 

When At = in equation (4), the secondary peak, 
Cs{p — (t>hi)^ contains all the required information to 
determine the altitude and C'^{h) strength for a given 
layer. It has been shown that Cs{p — (j^hi) is approxi- 
mately the difference between parallel and perpendic- 
ular slices from the 2D covariance with respect to the 
direction of the double star, Cb,\\ and Cb,± respec- 
tively, and can be written as (Avila et al. 1997) 



Cs{p - (l)hi) 



Cs,||(p,0,O)-Cb,±(p,0,O) (6) 



^0 



K{p, hi)CN{hi)dhi + n(p). 



where n{p) is the measurement noise, resulting in a sin- 
gle covariance peak shifted by (j)hi. (Note the change in 
subscripts.) K{p,hi) is the theoretical spatial covari- 
ance of the scintillation from a single star produced by 
a layer at height hi where f C^{hi)dhi = 1, assuming 
Kolmogorov turbulence. 

The altitude of a given later can be found from the 
spatial sampling across the aperture, Ar, and the stel- 
lar separation, 0, where the altitude resolution of the 
system. Ah, can be defined as (Kliickers et al. 1998) 



Ah- 



Ar 

) sec ( ' 



(7) 



C is the angle from zenith that the measurement was 
taken at. For Ar = 0.01 m.pix"^ and = 4 arcsec- 
onds, Ah ~ 500 m.pix"^ at zenith. 



Using the methods discussed in Johnston et al. 
(2000) and Johnston et al. (2002), C%{h) profiles can 
be obtained from UC-SCIDAR data. 

Using At > 0, it is assumed that if a layer moves 
with a horizontal wind velocity V{hi), then the scin- 
tillation pattern produced at the ground would move 
with the same V{hi). For a single turbulent layer, two 
scintillation patterns separated in time by At would 
be separated by V{hi)At (Avila et al. 2001). Hence 



CB{V{h^)At, (t), At) ^ Cs(0, 0, 0), 



(8) 



such that the height of the auto- covariance peak (i.e. 
At = 0) is approximately the height of the spatio- 
temporal cross-covariance peak displaced by a distance 
V{hi)At. To adequately capture motion of slowly mov- 
ing layers a long At is required. However, long At val- 
ues are blind to rapidly moving layers. Using multiple 
values of At, a full V{h) profile can be obtained. The 
algorithm used to analyse spatio-temporal covariances 
for V{h) profiles is presented in Mohr et al. (2008a). 



3 Data For Trending 

The first SCIDAR measurements taken at MJUO oc- 
curred in April 1999 using the system designed by Im- 
perial College (Johnston et al. 2002). These data re- 
vealed the presence of strong NGT and two different 
high altitude layers located at approximately 11 and 
13 km above sea level, with an estimated tq of 12.3 cm 
for the full profile. Temporal analysis indicated that 
the velocities of the 11 and 13 km layers were 6.45 
and 11.63 ms~^ respectively, whereas the NGT layer 
was attributed to dome seeing. However the measure- 
ments were collected over a single observation run of 
10 nights where only 50% of the nights provided useful 
observing conditions (Johnston 2000). For the devel- 
opment of AO for MJUO, it was decided that a more 
complete picture was required of variations in the tur- 
bulence profile with respect to season and weather. 

SCIDAR measurements for trending purposes were 
collected at MJUO from 2005 - 2007 using 
UC-SCIDAR. The majority of measurements were 
taken using the 1-m McLellan telescope at a focal ra- 
tio of F/13.5. Due to the noise characteristics of the 
CCD cameras only a handful of double or binary star 
systems were suitable. The location of MJUO per- 
mitted the binary star systems a Cru and a Cen to 
be used for a large portion of the year. However in 
the summer months, when both a Cru and a Cen are 
too far from zenith, fainter star systems such as Eri 
and V Car could be used with the latest version of 
UC-SCIDAR, which has more sensitive detectors. For 
the UC-SCIDAR system stellar separations should be 
limited to between 4 and 20 arcsec with the appar- 
ent magnitude of the primary star no fainter than 3.5, 
where the magnitude difference should be limited to 
r-^ 2.5. Table 1 gives the stellar parameters for the 
stars for measurements presented in this paper. 

Not all of the data collected between 2005 - 2007 
was suitable for use in the profiling of turbulence at 
MJUO. Only data with exposure times < 3.5 ms were 
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Table 1: Stellar parameters for stars used with 
UC-SCIDAR data presented. 



Star 


(j) (arcsec) 


mi 


Am 


a Cen 


8.7- 13.3 


-0.01 


1.36 


a Cru 


3.9 


1.25 


0.3 


Eri 


8.4 


3.20 


0.92 


V Car 


5.0 


3.02 


2.98 



included in the trending analysis, to ensure that turbu- 
lence was not subject to excessive blurring, resulting in 
the underestimation of layer strengths. Kliickers et al. 
(1998) used exposure times raging from 1.6 - 2.7 ms, 
whereas Avila et al. (2008) report use of exposures of 
3 ms. Measurements taken at zenith angles, greater 
than 35° were excluded from site profiling. Measure- 
ments with 30° < C ^ 35° were included only if sup- 
porting data taken at < 30° was acquired within 
approximately 30 minutes of the run. This was to en- 
sure that a suitable altitude range was sampled and 
that measured NGT was not subject to vertical air 
flows located near ground. A significant portion of the 
data collected during 2006 was subject to corrupt CCD 
readout and hence was not included in this study. 

Each observation period usually consisted of three 
to four consecutive nights. Data sequences typically 
consisted of 5000 frames from each camera which was 
recorded into file blocks of 500 frames per file. To 
decrease the processing time, 2D spatio-temporal co- 
variances for each camera were calculated using 2500 
frames with time delays between consecutive frames of 
At = 0, dt, 2dt . . . 6dt, where dt represents the frame 
rate of the CCD camera used. This ensured that at 
least 1000 cross-correlations were used in the longer 
At ensembles. It should be noted that all available 
frames were used such that for a At = 3dt, frame 1 
was correlated with frame 4, frame 2 with frame 5, 
frame 3 with frame 6, and so on. 

Measured velocities were cross-checked between se- 
quential runs and the various At covariances, and layer 
heights from temporal analysis were cross-checked with 
the corresponding Cj^{h) analysis for a given run to 
eliminate any falsely detected layers. Data collected 
in 2007 was collected at 60 Hz, whereas data collected 
prior to 2007 was collected at 30 Hz due to a limitation 
in the CCD cameras used at the time. 

In a significant number of cases there was simulta- 
neous pupil-plane and generalised SCIDAR data. This 
provided an added check to layer heights in the free at- 
mosphere, as well as an insight into how the strong low 
altitude layers can affect the measurements obtained 
for high altitude layers. 

Across the three years, a total of 324 pupil-plane 
and 315 generalised SCIDAR data sequences were suit- 
able for C%{h) profiling. 134 pupil-plane and 112 gen- 
eralised SCIDAR data sequences were used for V(h) 
profiling. The significant decrease in the number of 
suitable runs for V{h) profiling was due to the frame 
rate of early UC-SCIDAR data and the nature of NGT 
present at the time in later UC-SCIDAR data. Strong 



NGT can blur any covariance peaks detected during 
temporal analysis, particularly when using star sys- 
tems with a narrow angular separation. As such it is 
not always possible to determine the height of layers 
using temporal analysis as the primary and secondary 
peaks blur together. 

The C%{h) profiles presented in this paper include 
dome/mirror seeing as any AO system developed for 
the McLellan 1-m telescope at MJUG would need to 
compensate for dome/mirror seeing too. 

All ro, Oo and fc values calculated from 
UC-SCIDAR data have been determined for a wave- 
length, A, of 589 nm, although the measurements were 
collected using broadband white light. The variable 
weather at MJUO (ranging from calm, clear nights 
to gusting winds and thickening clouds) resulted in 
a variety of profiles being detected. During increas- 
ing cloud cover, with moderate to high ground wind 
speeds, the NGT present was exceptionally strong re- 
sulting in a pupil-plane ro that was similar to that 
found for the corresponding generalised data during 
analysis. In these cases the generalised ro estimate 
was used in the determination of the pupil-plane av- 
erages to remove any bias toward noise in pupil-plane 
results. Results are presented from the most recent to 
the earliest. 

4 Trends at MJUO 
4.1 Cjf{h) Profiles 

Figures 4 and 5 show the pupil-plane and generalised 
C]^(h) profiles respectively for data collected in 2007. 
The NGT measurements at MJUO tend to dominate 
and typically mask any activity present in the upper 
layers, as seen in Figures 4(a) and 5(a). To reveal 
possible features in these upper layers the C%{h) pro- 
files are scaled so that the colour range is limited to 
C%{h)Ah values between 10"^^ and 10"^^ m^/^ Val- 
ues below 10~^^ m^/^ are likely to result in images that 
are diffraction-limited rather than turbulence-limited 
on a 1-m telescope, ro is 1.55 m for a layer with 
C%{h)Ah = 10"^"^ m^/^ and a wavelength of 589 nm. 
Values greater than 10~^^ m^^^ are approaching lev- 
els that are classified as strong turbulence (Andrews 
2004). Figure 4(b) shows the colour scaled images 
for pupil-plane measurements from 2007. Figure 5(b) 
shows the colour scaled profiles for 2007 generalised 
measurements. 

Examination of the overall pupil-plane C%{h) mea- 
surements (Figure 4) indicates that there was signifi- 
cant turbulence located at low altitudes (i.e. below 5 
km above sea level) with an additional layer that could 
be seen in some measurements at 12 - 14 km above sea 
level. Where low altitude turbulence was strong (i.e. 
all of January, the latter half of May and the major- 
ity of June) the upper altitude layers were masked. In 
the generalised data (Figure 5(b)) little to no upper 
altitude activity was detected. The only exception oc- 
curred in the first half of May where little NGT was 
detected in the pupil-plane data. However the strength 
of the detected layer at approximately 14 km above sea 
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(b) Scaled Pupil-plane 



Figure 4: Pupil-plane C'^{h) profile trends ob- 
served over 2007. Gaps have been added where no 
data is present for more than two hours. Data from 
the individual months is separated by solid black 
lines. The image in (b) has been scaled such that 
any C'^{h)Ah value above 10~^^m^/^ is set to the 
maximum colour range and any C'^{h)Ah value 
below 10~^^m^/^ is set to the minimum colour 
range. 



level does not match the strength of the layer detected 
at comparable height in pupil-plane data. 

From equations (4) and (6), there is an assump- 
tion that each layer of turbulence present above a site 
is statistically independent and hence the strength of 
the layers found is not dependent on other layers in 
the structure. If this assumption holds true then the 
measured strength of the detected turbulence from any 
given high altitude layer should be the same regard- 
less of whether pupil-plane or generalised measure- 
ments were employed. However issues arise under con- 
ditions of medium to strong NGT resulting in an un- 
derestimate of the strength of turbulence located in 
the higher levels in generalised SCIDAR measurements 
(Mohr 2009). A correction factor can be applied to 
the high altitude layers detected in generalised data. 

Figure 5(c) shows the colour scaled corrected gen- 
eralised data for 2007. The correction factor used was 
a localised Gaussian curve centred around the alti- 
tudes of the layers detected in the pupil-plane data, 
correcting for differences in measurement plane, with 
a standard deviation of 2 Ah, where Ah is the altitude 
resolution of the measurement. The level of correction 
needed is dependent on the altitude of the given layer 
and the strength of the NGT detected, based on sim- 



. 15 

CD 
CO 

CD 

g 10 
< 

CD 

3 5 

s 





h 

li s ^ 



2007 Jan 2007 May 

Month of Observation 



(a) Generalised 



























hi 1 1 i 


A, II 




2007 Jan 2007 May 2007 Jun 
Month of Observation 

(b) Scaled Generalised 


























hi \ M i 


ii. 1 i 





2007 Jan 2007 May 

Month of Observation 



(c) Scaled Corrected Generalised 

Figure 5: Generalised C'^{h) profile trends ob- 
served over 2007. Time and C'^{h) colour scal- 
ing used is as per Figure 4. (c) Correction factors 
applied to turbulent layers detected in the free at- 
mosphere is incorporated. See text for details. 
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Figure 6: Pupil-plane and generalised C'^{h) pro- 
file trends observed over 2005. Time and C'^{h) 
scaling, and generalised data corrections used are 
as per Figure 5. 



Figure 7: Pupil-plane and generalised C'^{h) pro- 
file trends observed during autumn from 2005 to 
2007. Time and C%{h) scaling, and generalised 
data corrections used are as per Figure 5. 



ulations (Mohr 2009). It should be noted that this 
correction has little effect on the estimated ro for the 
profiles, due to the strength of the NGT layer that 
dominates the profiles. Following the correction it can 
be noted that similar strength turbulent layers are seen 
in both the pupil-plane and generalised SCIDAR data 
for the high altitudes (Figures 4(b) and 5(c)). Where 
possible a similar correction has been applied to all 
generalised data presented in this paper. 

The layer found at 12 - 14 km above sea level can 
be associated with turbulence found in the tropopause 
region, which is commonly incorporated into models 
for C%{h) profiles (Hardy 1998). For a site such as 
MJUO, a significant level of turbulence at low alti- 
tudes is expected. MJUO is located roughly 50 km 
west of the Southern Alps. With prevailing westerlies 
over much of New Zealand, the low- to mid- altitude 
wind structure is significantly affected by the terrain 
associated with the Southern Alps (Sturman and Tap- 
per 1996). 

For data collected in 2005, shown in Figure 6, a 
strong low altitude layer located at less than 5 km 
above sea level is visible. Also present is a weaker high 
altitude layer that ranges between 10-14 km above 
sea level. In some months an additional layer can be 
seen at 6 - 8 km above sea level. 

Figure 7 shows the scaled pupil-plane and gener- 
alised Clj{h) profiles obtained for the autumn months 
(i.e. April and May) for UC-SCIDAR data collected 
in both 2005 and 2007. During 2005 April and May, a 
weak high altitude layer was found at approximately 



11 - 12 km above sea level with a C%{h)/S.h strength 
of approximately 3 x 10"^^ m^/^. In 2007 May the 
height of this high altitude layer increased by approx- 
imately 2 km. The height of the tropopause is known 
to fluctuate up to 4 km throughout the year (Sturman 
and Tapper 1996). It is reasonable to assume that 
this is turbulence generated at the same pressure scale 
height. NGT and boundary layer turbulence is seen 
to extend up to 5 km above sea level. This is seen in 
both the pupil-plane and generalised data. Little can 
be ascertained from the generalised data as to a trend 
in the height and strength of the high altitude layer, 
due to the strength of the NGT and the noise present 
in the 2005 data. 

Figure 8 shows the scaled C]^(h) profiles obtained 
for the winter months (i.e. June and July) with UC- 
SCIDAR collected in both 2005 and 2007. A high al- 
titude layer is seen consistently at approximately 12 
km above sea level with an average C%{h)Ah strength 
of approximately 3 x 10~^^ m^^^. This is of similar 
strength to that seen in the autumn months (i.e. April 
and May) , suggesting little to no difference in the high 
altitude layer. However, an additional mid-altitude 
layer is seen in some data, with heights ranging from 6 
- 8 km above sea level and with varying strength. As 
seen with the April/May data, low altitude turbulence 
extends up to 5 km above sea level. Again little can be 
ascertained from the generalised data about trends in 
the heights and strengths of the high altitude layers. 

The monthly averages and standard deviations for 
ro and Oq are summarised in Table 2. Also shown is 
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Figure 8: Pupil-plane and generalised C'^{h) pro- 
file trends observed during winter from 2005 to 
2007. Time and C%{h) scaling, and generalised 
data corrections used are as per Figure 5. 



the mean turbulence altitude, /lo, which can be found 
using (Garcia- Lorenzo et al. 2009) 



/lo = 0.314 



(9) 



providing a measure of the effective height of dom- 
inant turbulence if the multi-layer structure was re- 
placed with a single layer. Note that /lo values shown 
reflect the distances above sea level. For generalised 
data, ho has been corrected for defocus distances for 
each observation sequence. 

January 2007 data exhibited similar ro and val- 
ues for both pupil-plane and generalised data, which 
indicates a significant low-altitude layer. This is also 
shown in the Hq values of 4.8 and 2.9 km for pupil- 
plane and generalised data respectively. In summer, 
the longer, warmer days heat the surrounding ground 
and buildings which can lead to greater NGT effects. 
The pupil-plane Oq values obtained for May are smaller 
than those obtained for June although the ro values 
are similar for the two months. The difference in Oq 
can be attributed to the high altitude layers found in 
May, which were stronger and higher than those found 
in June (Figure 5(b)). This is also seen in the higher 
pupil-plane Hq value for 2007 May of 7.3 km compared 
to 5.7 km for 2007 June. 

During 2005, the winter months (i.e. June and 
July) have a pupil-plane fo in the order of 20 cm, 
whereas autumn (i.e. March, April and May) and 
spring (i.e. August) have a pupil-plane fo ranging from 
approximately 10 cm to 15 cm. ^o follows a similar 



trend to that of fo, with the winter months having 
a value of approximately 2 arcsec compared to 1.3 - 
1.7 arcsec for autumn and spring, ^o for 2005 March 
was 1.6 arcsec with a large standard deviation gqq of 
0.7. This is associated with the noise present in a sig- 
nificant portion of the runs during 2005 March. The 
fo values for generalised measurements are reasonably 
consistent throughout the two years; 5-6 cm for 2007, 
and 7 cm for 2005. Little variation is seen for gener- 
alised ^0 values across the entire campaign, suggesting 
a consistent, dominating NGT layer. The significant 
variations seen in pupil-plane fo and Oq values suggest 
that the strengths of the high altitude layers fluctuate. 
The average ro for all UC-SCIDAR data was 12 ± 5 
cm and 7 =b 1 cm for the pupil-plane and generalised 
measurements respectively. Oq was 1.5 ±0.5 arcsec and 
l.OzbO.l arcsec for the pupil-plane and generalised mea- 
surements respectively. Hq was 6 ± 1 km and 2.0 ± 0.7 
km for pupil-plane and generalised measurements re- 
spectively. The values of Hq obtained does reinforce 
that the dominating turbulence at MJUO is located 
at near ground altitudes. 

UC-SCIDAR estimates for ro are smaller than those 
obtained in 1999 April using the Imperial College sys- 
tem (i.e. 12.3 ±1.4 cm) (Johnston et al. 2002). Al- 
though layer height estimates were similar for the two 
systems, the strength of the NGT layer was in the or- 
der of 6 - 10 times stronger for UC-SCIDAR data. 
Some variation can be expected due to the amount of 
time that has passed between the two instruments be- 
ing used. However the lower ro values obtained using 
UC-SCIDAR do match the observation conditions typ- 
ically seen by observers, which has a nominal angular 
resolution, ^res, of ^ 2 arcsec (A. Gilmore (MJUO) 
2006, private communication), ^res from UC-SCIDAR 
was 2.5 arcsec for the full profile, calculated at a wave- 
length of 589 nm. ^res for data taken in 1999 April was 
1.2 arcsec. 

The large variation in the pupil-plane measure- 
ments suggests not only a possible relationship with 
seasonal changes, but also with the weather at the site. 
June and July of 2005 saw high pupil-plane fo and Oq 
values which can be attributed to the calmer weather 
seen during these observational periods. However most 
other months show similar fo and Oq estimates within 
the margin of error. This suggests that weather con- 
ditions have a greater influence on the profiles ob- 
tained. It is suggested that data obtained from UC- 
SCIDAR be correlated to meteorological data to inves- 
tigate whether SCIDAR data could predict weather- 
related seeing over the site. 

4.2 V{h) Profiles 

Figure 9 shows the average wind speed, |V^(/i)|, mea- 
sured as a function of height, /i, for the 2007 data. 
There is some scatter, due to errors associated with 
height estimation, but the measurements show a tur- 
bulent layer at 11 - 14 km above sea level with an 
average speed of 18 ms~^, with velocities ranging from 
^ 6.5 ms~^ to over 30 ms~^. In some cases a mid- 
altitude layer at 6 - 7 km above sea level can be seen 
moving at approximately 7 ms~^. On nights with sig- 
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Table 2: Monthly averages and standard deviations for ro, and ho for all months, ho values are indicated 
as distances above the sea level. Elevation of MJUO is 1024 m. Generalised values have been computed 
based on corrected profiles. Generalised values of ho have been corrected for defocus distances. 



Month 


^0 

(cm) 


Pupil-plane 

O'ro Oo 

(cm) (arcsec) 


(arcsec) 


ho 
(km) 


^0 

(cm) 


Generalised 

(7rQ Oo 

(cm) (arcsec) 


(arcsec) 


ho 
(km) 


2007 June 


12 


3 


1.8 


0.7 


5.7 


5 





1.1 


0.3 


2.3 


2007 May 


11 


7 


1.3 


0.9 


7.3 


6 


1 


1.0 


0.2 


2.6 


2007 January 


6 


3 


0.8 


0.2 


4.8 


5 


2 


0.9 


0.1 


2.9 


2005 August 


10 


2 


1.7 


0.4 


4.5 


8 


1 


1.1 


0.1 


1.6 


2005 July 


18 


4 


2.2 


0.4 


6.3 


7 


1 


1.2 


0.2 


1.3 


2005 June 


22 


3 


1.8 


0.3 


8.8 


6 


1 


1.0 


0.1 


1.5 


2005 May 


12 


1 


1.7 


0.2 


5.0 


7 


1 


0.9 


0.2 


1.1 


2005 April 


11 


3 


1.3 


0.3 


6.8 


7 


1 


0.9 


0.1 


1.9 


2005 March 


15 


3 


1.6 


0.7 


7.2 


8 


2 


1.1 


0.3 


2.8 


1999 April 












12.3 
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Figure 9: Average wind speeds, |1^(/?')|, for obser- 
vations taken during 2007. Measurements from 
the individual months are separated by solid black 
lines. 



nificant NGT, a low-altitude layer was seen with wind 
speeds ranging between 10 to 24 ms~^. 

Velocity measurements obtained in April 1999 in- 
dicated high-altitude layers at 11 and 13 km travelling 
at 6.45 and 11.63 ms~^ respectively. This is consistent 
with the measurements obtained using UC-SCIDAR 
for layers at this altitude. 

Table 3 shows the monthly averages and standard 
deviations for fc- Note that the value shown for 2007 
January is averaged across both the pupil-plane and 
generalised data. The values for 2007 June reflect the 
average for data collected on May 31 and June 1 only 

The monthly calculated Greenwood frequency av- 
erages, /g, for 2007 were approximately 30 Hz or less. 
For a wavelength of 589 nm, the estimate for fo ranges 
between 30 - 90 Hz depending on the models used for 
C%{h) and V{h) profiles (Tyson & Frazier 2004). The 
measured fc obtained at MJUO is at the bottom end 
of this range. It is most likely that the fc has been 
underestimated due to the gaps that exist in the V{h) 
profiles. 

fc calculated from SCIDAR can be underestimated 
because wind velocity measurements are reliant on co- 



Table 3: Nightly averages for fc for all months. 
Generalised values have been computed based on 
corrected profiles. 





Pupil- 


-plane 


Generalised 


Month 






fc 






(Hz) 


(Hz) 


(Hz) 


(Hz) 


2007 June 


26 


12 


29 


31 


2007 May 


30 


23 


8 


13 


2007 January 


20 


15 






2005 July 


10 


7 


11 


14 


2005 June 


13 


5 


10 


18 



variance strengths in the measurement plane being suf- 
ficiently strong with respect to the background covari- 
ance noise. Aperture normalisation in the data am- 
plifies noise, especially near the aperture edge. Peaks 
approaching the aperture edge can be hidden by the 
noise. The detection of only the central peak with one 
secondary peak, termed partial triplet analysis (Mohr 
et al. 2008a), does allow for more layer velocities to 
be found. In addition, some covariance peaks may be 
obscured due to their close proximity to other peaks. 
Although all layers seen in V{h) will have an associ- 
ated Cj^{h) strength, not all C^(h) layers will have a 
measurable V{h) due to the position of the covariance 
peaks relative to other peaks and the aperture edge 
and the resulting covariance strength. This results in 
gaps in the measured V{h) profile and is reflected in 
the large af^ values. 

The UC-SCIDAR system used in 2005 captured 
data at a frame rate of 30 Hz, which limits the maxi- 
mum detectable velocity on a 1-m telescope to 30 ms~^ 
under ideal conditions (Mohr et al. 2008b). This was 
a limitation of the CCD cameras used at the time. 
Based on the V(h) profiles from the current system, 
which has a frame rate of 60 Hz, data collected dur- 
ing 2005 should provide a reasonable V{h) profile for 
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Figure 10: As per Figure 9, except for 2005 data. 
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Figure 11: As per Figure 9, except for observations 
taken during winter months of 2005 and 2007. 



low- to mid-altitude layers, however high-altitude layer 
velocities may not be measurable. Temporal analysis 
was performed on data from June and July 2005 data 
only. These months provided the most reliable data 
set with a large number of runs utilising exposures of 
1 - 2 ms. 

Figure 10 shows the average wind speeds for 2005 
June and July. In June a layer at approximately 12 
km above sea level was moving consistently at 12 - 
15 ms~^. In July there was much more scatter in the 
average velocity. 

Using the V{h) profiles obtained for 2005 June and 
July, the average Greenwood frequency, fo , was found 
to range between 10 and 20 Hz. There are significant 
gaps present in the V{h) profiles, as indicated by the 
standard deviation in /g, ct/^, ranging between 5 and 
15 Hz (Table 3). As such fo is likely to be underesti- 
mated. 

Figure 11 shows the average wind speed for the 
winter months (June and July). Although a significant 
amount of scatter exists for the data from 2005 July 
and 2007 June the layer heights found are similar. The 
high altitude layer has an average speed of 10 - 15 
ms~^ in calmer weather, but speeds in excess of 25 
ms""*^ at other times. 



5 Finding a Model that Fits 

5.1 Cj^{h) Model 

Optical turbulence is highly irregular, where turbu- 
lence strengths can vary by an order of magnitude 



around a mean value (Hardy 1998). Commonly used 
models represent a mean profile of C^(/i) measure- 
ments taken over extended periods of time. 

The Hufnagel- Valley (HV) model is commonly used 
to describe the average turbulence at an astronomical 
site (Tyson & Frazier 2004). A standard HV model 
consists of three main components: an exponentially 
decreasing C^(/i) through the troposphere; a peak at 
approximately 10 km above ground corresponding to 
a tropo-pause layer; and a strong surface layer (Hardy 
1998). 

At many sites additional layers have been detected 
at low- to mid-troposphere altitudes (Avila et al. 2004; 
Fuensalida et al. 2004; Prieur et al. 2001; Wang et al. 
2008). The generic HV model, incorporating an addi- 
tional mid- altitude layer, is a sum of exponential terms 
such that C%{h) is given by (Hardy 1998) 



Cl{h) 



A exp 
-\-B exp 



-) 

Ha J 

h_ 

Hb 



+^exp(-(^), (10) 

where A is the turbulence coefficient for near-ground 
turbulence (i.e. oc C]^{0)) and Ha is the height for its 
1/e decay, B and Hb are similarly defined for turbu- 
lence in the troposphere, and C and He are related to 
the turbulence peak located at the tropopause. The 
fourth term in equation (10) can be used to define one 
or more isolated layers, where D and Hd define the 
strength and height of the layer and d specifies the 
layer thickness. 

Figure 12 shows selected pupil-plane and gener- 
alised (h) profiles obtained from various months in 
2005 and 2007. Also shown are three different models 
all based on the HV model. Table 4 lists the param- 
eters for the various models shown in Figure 12. The 
HV 5-7 model (indicated by the solid black line) has 
parameters such that the resulting ro and Oq are 5 cm 
and 7 //radians (1.44 arcsec) when using A = 500 nm. 
The HV 5-7 model, while producing ro and Oq values 
appropriate for the site (i.e. 6 cm and 1.7 arcsec re- 
spectively for A = 589 nm), results in a tropopause 
layer that is slightly too low, residing at 10 km above 
the site (i.e. 11 km above sea level) and too weak. By 
increasing the altitude of the tropopause peak by 500 
m and increasing the coefficient C from 3.59 x 10~^^ 
to 5.94 X 10~^^ (MJUOl model) then the Oq decreases 
to 0.96 arcseconds which is more in line with the mea- 
sured values from the generalised C^{h) profiles. 

The near-ground and low-altitude turbulence is seen 
to regularly extend up to 5 km above sea level (Fig- 
ures 4-8). An isolated layer was added with a peak 
at 1.5 km above the site (i.e. 2.5 km above sea level) 
and a thickness of 1 km. The peak strength of this 
layer has been set to 2 x 10~^^ m~^/^. The modified 
HV model, with the addition of this low-altitude layer, 
is called MJU02. In some profiles an additional layer 
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was found at heights ranging from 6-8 km above sea 
level. MJU03 incorporates an additional layer at 6.5 
km above sea level, and is the recommended model for 
use in the AO design for MJUO. The profiles of the HV 
5-7 model, MJUOl and MJU03 are shown in Figure 
13. 

Additional layers at approximately 2 and 6.5 km 
about sea level have also be regularly detected at other 
sites. Avila et al. (2001) reported the detection of five 
different layers during the 1998 site testing campaign of 
Cerro Tololo and Cerro Pachon in Chile. Strong layers 
were found at the tropopause region (i.e. 11-12 km) 
and at low altitudes, extending up to approximately 4 
km above sea level. Although the layer found at 6.5 km 
was often weaker than the low altitude turbulence, it 
was consistently present. Similar low- to mid- altitude 
turbulence was detected at San Pedro Martir, Mexico 
(Avila et al. 2007). 

The HV model is such that the resulting peak 
present at the tropopause region is not a spike but 
rather a broad peak, which accounts for variations in 
layer heights and fluctuations in the turbulence strength 
seen over time. Instantaneous spikes and variations in 
the heights of the turbulent layers are smoothed out. 
This has the effect of broadening the peaks seen, par- 
ticularly in the tropopause region. Although it would 
be ideal to refine the model such that the breath of 
the tropopause region is not so wide, as a model for 
AO design MJU03 is the recommended model. Using 
a wavelength of 589 nm, tq is estimated to be 6 cm for 
MJU03. Oo is estimated to be 0.9 arcseconds. ho is 
estimated at 4.2 km above the measurement plane. If 
a defocus distance of 3 km was used the resulting ho of 
1.2 km would be in line with the values obtained from 
the generalised UC-SCIDAR measurements. 




1 2 3 4 5 6 

C^(h)Ah xlO"^' 

(a) Selected pupil-plane and generalised data 




10" 



10 



10 10 

C^(h)Ah (m"^) 

(b) Semi-logarithmic plot of (a) 



Figure 12: C'^{h) model fitting. Both (a) and 
(b) show selected pupil-plane and generalised data 
from (*) 2007 June, ( _ ) 2007 May, (□) 2005 June 
and (0) 2005 April. Models shown are (— ) the HV 
5-7 model, (- -) MJUOl: a modified HV model 
and (- • -) MJU03: a modified HV model that 
incorporates two additional layers. 



5.2 V{h) Models 



To describe the wind velocity with increasing altitude, 
a Greenwood wind model is commonly used (Hardy 
1998). The Greenwood model is a Gaussian based 
model and is defined as (Tyson & Frazier 2004) 



V{h) - U(0) + V{Ht) exp 



h cos ( — Ht 



Lt 



[sin^ P + cos^ P cos^ 



2 >ll/2 



(11) 



where U(0) is the wind velocity at ground level, V{Ht) 
is the velocity at the tropopause located at an altitude 
i^T, Lt is the thickness of the tropopause layer, /3 
is the wind direction relative to the telescope azimuth 
and C is the angle of the telescope from zenith. The di- 
rection corrections used in the Greenwood wind model 
are for a telescope that uses a horizontal coordinate 
system and strictly speaking should not be applied to 
telescopes at MJUO as they employ an equatorial co- 
ordinate system. However for AO design it is the wind 
speed that is important. For the purpose of further 
discussion it will be assumed that /3 = 0° and hence 




Figure 13: C'^{h) turbulence models for MJUO. 
(— ) HV 5-7 model; (- -) MJUOl: Modified HV 
model; (- • -) MJU03: Modified HV model in- 
corporating two additional layers. MJU03 is the 
recommended model. 
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Table 4: Parameters for C'^{h) turbulence models using a generic HV model. 



Moder 


C 


He 


D 


Hd 


d 


to' 


Oo 




(xlO-^^) 


(m) 


(xlO-^^) 


(m) 


(m) 


(cm) 


(arcsec) 


HV 5-7 


3.59 


1000 









6 (27) 


1.74 (1.81) 


MJUOl 


5.94 


1050 









6(16) 


0.94 (0.96) 


MJU02" 


5.94 


1050 


2 


1500 


1000 


5 (16) 


0.92 (0.95) 


MJU03^ 


5.94 


1050 


2 (0.3) 


1500 (5500) 


1000 (500) 


5 (15) 


0.90 (0.94) 



^For all models indicated A = 17 x IQ-^^, Ha = 100 m, B = 27 x 10'^^, and Hb = 1500 m. 

^ro and values are specified for A = 589 nm for the full profile and for /i > 3 km above the telescope in brackets. 
^MJU02 includes a strong low-altitude layer. 

^MJU03 incorporates an additional mid-altitude layer for which the parameters are indicated in brackets. 



equation (11) becomes 



V{h) = y(0)+y(i/T)exp 



h cos C — Ht 



xcos C- 
(12) 

Figure 14 shows the V{h) measurements from 2007 
May 3. The V{h) profile for these data has an ideal 
profile that can be modelled by a modified Bufton wind 
model, where the Bufton model is a specific Green- 
wood model where V{0) — 5 ms~^ and V{Ht) — 30 
ms"^ with Ht 9.4 km and Lt 4.8 km for C = 0°. 
Also shown in Figure 14 is the Bufton and modified 
Bufton wind models. The standard Bufton wind model 
(indicated by the dashed green line) assumes that the 
tropopause layer is at 9.4 km above the telescope with 
a wind speed of 30 ms~^. This is low for the V{h) pro- 
file for May 3. Moving the model tropopause height 
to 12 km above sea level (indicated by the solid black 
line) allows the model to encompass the lower turbu- 
lent layers detected, as well as the activity detected 
in the tropopause region. If a zenith angle was in- 
corporated into the model (red dashed line) then the 
effective wind speed in the tropopause region decreases 
to a value that is better suited to the data. The modi- 
fied Bufton model used incorporated a zenith angle of 
20°, as this was the angle used by most of the mea- 
surements made on May 3. 

Regardless of the model parameters used, a Gaus- 
sian based model will not adequately describe the V{h) 
profiles detected when low altitude layers at 2 - 5 km 
have wind velocities between 15 - 20 ms~^. A large 
portion of the data collected for velocity had the pres- 
ence of some low-altitude wind. It was decided to add 
a second Gaussian peak, such that 



v[h) = y(o) 

+y(i7T)exp 
+y(ifi)exp 



h cos C — Ht 
Lt 

h cos ( — Hi 



(13) 
X cosC 

X cosC, 



where V{Hi) is the velocity of a low- altitude layer lo- 
cated at Hi above the telescope with a thickness of 
Li. 

Unlike the modelling of C]^(h)^ it is difficult to 
employ a generic model that would encompass the ma- 
jority of conditions, as the velocity seen in the upper 




10 20 30 

Wind Speed (m/s) 



I "° |V(h)|<10m/s □ 10<|V(h)|<20m/s □ |V(h)| > 20 m/s| 

Figure 14: Wind speed analysis for observations 
made on 2007 May 3. Models shown are ( ) the 
standard Bufton model, ( — ) a modified Bufton 
model with Ht = 11 km, and (- -) a modified 
Bufton model with Ht = 11 km and ^ = 20°. 
Elevation for MJUO is 1024 m above sea level. 



layers is dependent on the velocity detected near the 
ground. Figure 15(a) shows the instantaneous wind 
speeds obtained during June and July of 2005 and May 
and June of 2007. Overlaid are four different V{h) 
models developed to encompass the range of velocity 
characteristics detected. The parameters for the four 
models for MJUO are listed in Table 5. 

MJUO IV (indicated by the solid purple line in Fig- 
ures 15(a) and 15(b)) is designed for very calm and 
clear nights where little NGT is present and the high 
altitude turbulence has low wind speeds. MJU02V 
(indicated by dashed black line) was developed for 
nights such as 2007 May 3, where high altitude tur- 
bulence is strong, but low-altitude layer wind speeds 
are very low. Both MJUO IV and MJU02V are based 
on the traditional Greenwood wind model. MJU03V 
(indicated by the solid red line) and MJU04V (indi- 
cated by the dashed blue line) employ a second Gaus- 
sian peak described in equation (13). MJU03V is in- 
tended for moderate ground wind speeds measuring 
2.8 - 5.6 ms~^. Based on UC-SCIDAR measurements 
obtained this is the most likely situation to be encoun- 
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Wind Speed (m/s) 



(a) Fit of V{h) models to UC-SCIDAR data 




10 20 30 40 

Average Wind Speed (m/s) 



(b) MJUO V(h) models 

Figure 15: Fit of V{h) models to measured pro- 
files. Models shown are (— ) MJUO IV, (- -) 
MJU02V, (• • • ) MJU03V and ( ) MJU04V. 



tered and hence should be the preferred model for AO 
design. MJU04V was developed for situations where 
high ground wind speeds are present. 

Using the MJU03 Cj^{h) model and the MJU03V 
V{h) model for moderate ground wind speeds, fc is es- 
timated to be 79 Hz for a wavelength of 589 nm, which 
is more in line with the 30 - 90 Hz range suggested by 
Tyson & Frazier (2004). 

Although adequate models for V{h) profiles were 
determined, it is recommended that more temporal 
data be collected to refine the models. An investiga- 
tion into the correlation between low and high altitude 
velocities could also be conducted. 

6 Conclusions 

UC-SCIDAR measurements taken between 2005 - 2007 
detected strong NGT with a weaker layer located at 12 
- 14 km above sea level. On calm nights a third mid- 
altitude layer was detected at ^ 6 km above sea level. 
In a significant amount of data strong low altitude tur- 
bulence extended up to 5 km above sea level. The 
measurements suggest an average ro of 12 it 5 cm and 
7=bl cm for pupil-plane and generalised profiles respec- 
tively. This corresponds to an angular resolution, ^res, 
of 2.1 arcsec for the full profile. The average Oo values 
were 1.5 ±0.5 arcsec and 1.0 ±0.1 arcsec for pupil-plane 
and generalised profiles respectively. Average ho val- 
ues were 6± 1 km and 2.0 ±0.7 km for pupil-plane and 
generalised measurements respectively. 

During spatial analysis, correction factors were ap- 



plied to generalised profiles to the layers detected in 
the free atmosphere in the post-processing phase. It 
would be preferable to have these corrections incorpo- 
rated into the inversion algorithm to reduce the amount 
of double handling of data. To create an automatic de- 
tection system for this, a thorough investigation as to 
the effects of NGT strength on the detected strength 
from higher layers should be conducted, including var- 
ious altitudes and high altitude turbulence strength. 

Temporal analysis detected layers located at sim- 
ilar altitudes, with tropopause layer velocities of 12 - 
30 ms~^, dependent on weather conditions. Low al- 
titude turbulence layers had velocities ranging from 2 
ms~^ to well over 24 ms~^. No trends could be estab- 
lished for the values of /g, due to the gaps in the V{h) 
profiles. Little seasonal variation was detected in the 
C]^(h) profiles. Both Cj^{h) and V(h) profiles were 
highly dependent on the weather conditions. 

A modified Hufnagel- Valley (HV) model was de- 
veloped to describe the profiles, incorporating 
a strong NGT layer, a layer at 11 km above the tele- 
scope (i.e. 12 km above sea level) and two additional 
layers: one at 5.5 km above the telescope and the other 
at 1.5 km above the telescope extending up to 4 km. 
The resulting model estimates an ro of 6 cm for a wave- 
length of 589 nm, which corresponds to a ^res of 2.5 
arcsec. Oo is estimated at 0.9 arcsec. 

A series of V{h) models were developed, based on 
the Greenwood wind model with an additional Gaus- 
sian peak located at low altitudes to model the V{h) 
profiles seen at MJUO. The models correspond to calm, 
light, moderate and strong ground wind speed condi- 
tions seen at the site. Using the modified HV model for 
C%{h) profiles and the suggested model for V{h) pro- 
files in the presence of moderate ground wind speeds, 
fc was estimated at 79 Hz for a wavelength of 589 nm. 
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